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PREFACE 

One  problem  that  often  arises  in  practice  is  that  of  finding  gen¬ 
eralized  minimum  cost  multicommodity  network  flows.  This  problem  has 
not  previously  been  treated  in  the  technical  literature.  Its  impor¬ 
tance  is  clearly  demonstrated  by  its  many  applications.  Among  these 
are  the  military  applications  of  interdiction,  Military  Airlift  Command 
(MAC)  scheduling,  Logair  ,  and  inter-  and  intra-theater  transportation. 

Previous  work  on  multicommodity  flows  has  treated  networks  in 
which  individual  arcs  have  flow  capacities.  This  paper  generalizes 
that  notion  to  include  capacities  on  linear  combinations  of  flow  on 
several  arcs,  thereby  greatly  increasing  the  scope  of  the  problems 
that  can  be  treated.  Some  of  the  multicommodity  situations  treated 
here  which  have  not  been  handled  previously  are  the  allocation  of 
directional  flow  in  undirected  arcs;  node  capacities;  multimode  net¬ 
works  with  limited  vehicle  supplies;  and  the  efficient  use  of  scarce 
resources.  In  addition,  this  paper  proposes  a  new  solution  technique 
which  should  be  more  efficient  than  ones  given  previously. 

While  this  paper  is  directed  towards  a  technical  audience,  the 
results  will  be  of  value  to  those  working  in  the  applied  areas  of 
interdiction,  MAC  scheduling,  Logair,  inter-  and  intra-theater  trans¬ 
portation,  and  logistics  system  analysis. 

*Long  term  contract  airlift  service  within  continental  U.S. 
for  movement  of  cargo  in  support  of  military  services'  logistics 
systems. 
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SUMMARY 


This  Memorandum  presents  an  algorithm  for  1'inding  a  generalised 
minimum  cost  multicommodity  network  flow.  Previous  work  on  multi¬ 
commodity  flows  has  treated  more  limited  situations  £n  which  flows 
on  individual  arcs  are  capacitated,  or  upper  hounded.  No  problems 
have  been  treated  in  which  limited  resources  are  shared  among  several 
arcs,  instead  of  only  one.  However,  many  actual  network  flow  situations 
require  that  the  latter  be  dealt  with.  Some  general  examples  of  this 
are  operation  of  a  transportation  system  with  one  or  more  vehicle  types 
in  limited  supply,  operation  of  a  network  under  a  limited  budget  con¬ 
straint,  or  operation  of  a  network  where  undirected  arcs  (which  must 
be  represented  by  two  directed  arcs)  are  capacitated.  Specific  military 

if 

examples  are  the  allocation  of  aircraft  to  routes  for  Logair  or  the 
Military  Airlift  Command  (MAC),  where  air  corridors  have  virtually 
unlimited  capacity,  but  the  quantity  of  aircraft  and  crews  are  a  li¬ 
miting  factor.  Another  military  use  is  for  interdiction  where  users 
of  a  lines-of-communication  network  are  limited  by  their  vehicle  supply, 
but  not  by  actual  network  throughput  capacity.  This  problem  is  over¬ 
come  in  this  paper  by  the  introduction  of  generalized  constraints. 
Specifically,  upper  bounds  are  placed  on  flows  through  nonnegative 
linear  combinations  of  arcs  rather  than  only  through  individual  arcs. 

This  study  is  structured  as  follows: 

°  Section  I  discusses  the  generalized  minimum  cost  multi- 
commodity  flow  problem, 

o  Section  II  formulates  it  in  terms  of  its  arc  subset- 
chain  incidence  matrix, 

* 

Long  term  contract  airlift  service  within  continental  U.S.  for 
movement  of  cargo  in  support  of  military  services'  logistics  systems. 
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I.  INTRODUCTION 

Problems  that  frequently  arise  in  practice  are  those  of  finding 
generalized  minimum  cost  multicommodity  network  flows.  These  have 
not  previously  been  treated  in  the  technical  literature. 

Multicommodity  flow  problems  in  which  individual  arcs  are  sub¬ 
ject  to  flow  (capacity)  constraints  have  been  dealt  with  quite  exten¬ 
sively.  Ford  and  Fulkerson  [l],  Jewell  [2],  Sakarovich  [3],  and 
Saigal  [4]  have  treated  the  problem  of  maximizing  the  sum  of  all 
commodity  flows  subject  to  capacity  constraints  on  the  arcs  and  Tomlin 
[5]  has  treated  the  problem  of  finding  minimum  cost  multicommodity  flows 
under  the  same  constraints.  However,  this  formulation  is  often  in¬ 
adequate.  In  many  actual  network  flow  situations,  the  physical  capac¬ 
ities  of  individual  arcs  are  not  limiting,  but  resources  shared  by  units 
of  flow  on  several  arcs  are.  For  example,  when  one  wishes  to  assign 
flights  to  air  routes,  the  capacities  of  the  air  corridors  themselves 
are  virtually  unlimited;  however,  the  number  of  flight  assignments 
that  could  be  made  would  be  constrained  by  the  number  of  planes  and 
the  number  of  pilots  available.  Examples  of  assigning  flights  to  routes 
are  Logair  [6],  MAC  scheduling  [7],  and  commercial  airline  scheduling 
[8-9].  In  military  interdiction  problems,  users  of  lines  of  communica¬ 
tion  networks  are  often  limited  by  their  supply  of  one  or  more  vehicle 
types  rather  than  their  physical  arc  capacities  [10-11], 

A  general  situation,  which  includes  some  of  the  cases  discussed 
above,  is  one  in  which  the  physical  means  of  transport  or  vehicles  are 
in  short  supply.  Each  mode  of  transportation  would  then  impose  a 
joint  capacity  constraint  on  all  arcs  which  used  it.  Another  general 
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situation  is  operations  under  a  fixed  budget,  which  must  be  allocated 
to  units  of  flow  on  all  arcs.  Still  another  is  when  physical  arcs  of 
a  transportation  system  may  be  used  for  flow  in  either  or  both  direc¬ 
tions.  Such  an  arc  would  need  to  be  represented  by  two  directed  arcs 
and  its  capacity  would  be  an  upper  bound  on  the  flow  of  these  two  arcs. 

The  formulation  treated  in  this  paper  handles  all  the  above 
problems  in  addition  to  those  treated  previously.  Specifically,  cer¬ 
tain  subsets  of  arcs  are  identified.  Rather  than  an  upper  bound  being 
assigned  to  the  flow  of  each  arc,  an  upper  bound  is  assigned  to  some 
positive  linear  combination  of  the  arc  flows  of  each  subset.  These 
types  of  bounds  may  be  considered  as  joint  capacity  constraints. 
Individual  arcs  may  belong  to  zero,  one,  or  several  subsets.  Given 
these  constraints,  and  unit  flow  costs  on  individual  arcs,  it  is  re¬ 
quired  to  find  a  minimum  cost  multicommodity  flow,  Fulkerson  [l3-l^ 
and  Busacker  and  Gowan  [15]  have  treated  single  commodity  minimum  cost 
flow  problems  involving  one  joint  capacity  constraint.  The  present 
paper  has  the  advantage  that  more  than  one  joint  capacity  constraint 
may  be  treated  in  multicommodity  networks.  It  is  surprising  that 
while  the  scope  of  problems  handled  by  this  formrlation  is  an  order-of- 
magnitude  greater  than  that  handled  by  the  individual  arc  capacity 
formulation,  some  of  the  algorithms  which  have  previously  been  pre¬ 
sented  to  solve  the  latter  do  not  require  extensive  modifications  to 
solve  the  former. 

The  remaining  sections  of  this  paper  first  present  an  arc  subset- 
chain  formulation  of  the  problem,  special  cases  of  joint  capacity  con¬ 
straints,  and  a  method  of  solution  based  on  an  extension  of  Fulkerson's 


method  for  multicommodity  flows  with  arc  capacities  [1].  Later  a 
new  method  of  solution  is  demonstrated  that  should  be  more  efficient 
Finally,  special  cases  of  the  objective  (cost)  function  are  shown. 


Consider  a  network  [N,  A1  of  nodes  i  ■  1,  . N  and  directed 
arcs  connecting  certain  pairs  of  nodes  for  which  flows  of  Q  different 
commodities  are  to  be  assigned.  Designate 

(i,j)  the  directed  arc  beginning  at  node  i  and  ending 
at  node  j , 


the  cost  of  a  unit  of  flow  on  arc  (i,j) 
the  source  for  commodity  q 

the  sink  (or  terminal)  node  for  commodity  q. 


The  topology  of  the  network  is  specified  by  its  arcs  am1  nodes. 


Each  commodity  is  identified  by  its  source  (the  node  where  it  originates) 


and  its  sink  (the  node  where  it  terminates). 


A  chain  is  a  sequence  of  distinct  nodes  ij,  . . . ,  i  such  that 

(i,  »  i,  .,)  is  an  arc  for  k  <  n.  If  i.  *  s  and  i  =  t  ,  the  chain  is 

'  k*  k+1  1  q  n  q 

said  to  be  one  for  commodity  q, 

A  cycle  is  defined  as  a  chain,  except  that  i^  ■  i^.  Alterna¬ 
tively,  a  chain  or  cycle  may  be  defined  in  terms  of  its  arcs,  (i^, 
i^+j).  A  chain  or  cycle  is  positive  or  negative  depending  on  whether 
Mie  sum  of  its  arc  costs  are  positive  or  negative  respectively. 

Given  a  network  containing  no  negative  cycles,  it  is  required  to 

find  chain  flows  for  all  commodities,  the  number  of  which  shall  be 

designated  by  the  symbol,  Q,  that  minimizes  total  network  cost,  z,  sub¬ 
ject  to  the  m  nonnegative  linear  constraints,  (1.2),  where 
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(1.1) 


(1.2) 


E  kii  yu 


A.  y  .  ^  b  r  =  1,  ,  m 

ij  r  J  r 


and 


is  the  upper  bound  or  capacity  for  the  linear  flow  con¬ 
straint  r 

/V 

is  the  nonnegative  contribution  of  a  unit  of  flow  on 
arc  (i,j)  to  flow  constraint  r 
is  the  flow  on  arc  (i,j) 

Expressions  (1.2)  ma>  be  thought  of  as  joint  capacii"  constraints 


ij 


yii 


on  subsets  of  arcs,  ...,  Sm<  Arc  (i,j)  is  in  subset  r  if  and  only 
r  r 

if  A.  ,  >0  where  A.  .  is  the  number  of  capacity  units  for  subset  r  con- 
ij  ij 

sumed  by  a  unit  of  flow  on  (i,j)  and  b^  is  the  total  joint  capacity  for 
subset  r.  Any  arc  may  belong  to  zero,  to  one,  or  to  many  subsets.  For 
the  special  case  where  each  individual  arc  has  a  capacity,  the  coefficient 
matrix  for  (1.2)  would  be  an  identity  matrix. 

This  generalized  multicommodity  flow  problem  may  be  formulated  as 
an  arc  subset-chain  linear  program.  Specifically,  if  C^,  ...,  are 
all  the  chains  of  the  network  for  the  various  commodities  and  x^  the 
flow  on  chain  C  ^  ,  then  the  program  may  be  formulated  as  follows: 


* 

r 

The  nonnegativity  restriction  A. .  may  be  relaxed  as  follows. 
Suppose  that,  for  each  r,  if  (i,j)  weife  to  be  assigned  a  cost  of  Ar .  , 
the  total  cost  on  all  cycles  would  be  nonnegative;  then  all  results^ 
which  follow  in  this  paper  would  still  hold. 


t 
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(2.1) 

(2.2) 


Find  min  z,  xt  i  0  such  that 

z«  ^ctxt 

t^l 

n 

£  art*t  +  *n+r  *  br  . . 

t«l 


where 

(3.1) 


E 

(i.jTcCj. 


k .  . 
ij 


(3.2) 


(l.J)cCt 


Thus  the  rows  or  constraints  of  the  linear  program  correspond  to  the 
joint  capacity  constraints  on  the  arc  subsets  while  the  columns  or 


variables  (other  than  the  slacks,  x  , ,  ...,  x  )  correspond  to  the 

n-Tl  rrrm 

network  commodity  chains.  The  constraint  coefficient,  a^,  repre¬ 
sents  the  total  amount  of  joint  capacity  on  subset  used  by  a  unit  of 
flow  on  chain  Ct>  and  ct>  of  course,  is  the  total  cost  of  a  unit  of 
flow  on  chain  Ct> 
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III.  APPLICATIONS  OF  JOINT  CAPACITY  CONSTRAINTS 

Some  of  the  special  applications  involving  networks  with  joint 
capacity  constraints  are  those  involving  decisions  as  to  the  direction 
of  traffic  flow  on  streets,  networks  with  node  capacities,  general  re¬ 
source  constraints,  and  multimode  networks  with  vehicle  limitations. 

Directional  Traffic  Flow 

In  many  actual  network  situations,  arcs  are  undirected  with 
capacities.  Thus,  one  may  be  faced  with  a  traffic  situation  where  one 
would  want  to  determine  which  streets  should  be  one-way,  which  should  be 
two-way  and,  for  the  two-way  streets,  the  number  of  lanes  there  should 
be  in  each  direction.  However,  these  arcs  need  to  be  represented  by 
oppositely  directed  one-way  arcs.  If  one  is  dealing  with  single  com¬ 
modity  networks  with  nonnegative  arc  costs,  one  is  assured  that  an 
optimal  flow  pattern  exists,  using  at  most  one  of  these  arcs,  and  one  may 
therefore  assign  to  each  of  these  directed  arcs  an  individual  capacity 
equal  to  that  of  the  original  undirected  one.  However,  when  more  than 
one  commodity  is  involved,  situations  may  arise  where  any  optimal  solu¬ 
tion  involves  flow  in  both  directions  of  an  undirected  arc.  For  this 
reason,  capacities  on  the  individual  arcs  are  insufficient  to  describe 
the  undirected  arc  capacity.  However,  it  can  be  handled  by  a  joint 
capacity  constraint.  Specifically,  if  (i,j)  is  an  undirected  arc  with 
capacity  b_,  then  the  joint  capacity  constraint  for  the  two  directed 

arcs  replacing  it  is  x. .  +  x.,  ^  b.  .. 

Ji  ij 

Node  Capacities 

Nodes  as  well  as  arcs  often  have  limited  capacities.  For  example, 
one  might  be  limited  in  an  inventory  situation  by  the  capacity  of  a 
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storage  facility;  or  in  a  traffic  network  the  capacity  of  an  inter¬ 
section  may  be  more  limiting  than  the  capacity  of  all  arcs  leading 
into  it.  These  situations  may  also  be  phrased  in  terms  of  joint  arc 
capacities.  For  a  source  node,  this  could  be  expressed  as  the  sum  of 
flows  on  arcs  leaving  the  source  cannot  exceed  the  source  capacity;  for 
a  sink  node,  the  sum  of  all  flows  on  arcs  entering  it  cannot  exceed  the 
sink  capacity;  and  for  any  other  node,  its  capacity  may  be  expressed  as 
an  upper  bound  on  the  sum  of  flows,  either  on  arcs  entering  it  or  leaving 
it.  Thus  a  capacity  of  b^  on  an  intermediate  node,  i,  may  be  expressed 
either  as 


)x,  .  <  b.  or 

L-j  ij  i 


<  b  . 


l 


) 


while  the  capacity  of  a  source  node  would  be  the  former  and  a  sink  node 
the  latter  of  these  two  expressions.  Node  capacities  for  single  com¬ 
modity  networks  have  previously  been  treated  in  [12], 


Resource  Constraints 


Another  situation  that  often  occurs  in  practice  is  operations 
under  resource  or  budget  constraints.  Specifically,  there  exists  li¬ 
mited  quantities  of  one  or  more,  resources,  such  as  men,  vehicles,  dol¬ 
lars  or  other  appropriate  items.  Suppose  each  unit  of  flow  on  arc 
(i,j)  consumes  units  of  a  particular  resource;  and  b  is  the  total 
supply  of  that  resource.  Then  this  could  be  expressed  by  the  joint 

capacity  constraint  on  all  arcs  I  A. ,x, .  s  b.  Such  resource  constraints 
K  ij  ij 

frequently  occur  in  aircraft  allocation  problems,  including  both  military 


-9- 


airlift  such  as  Logair  [6]  and  the  Military  Airlift  Command  (MAC) 
scheduling  [7]  and  commercial  airlift  [8-9].  Here,  the  capacities  of 
the  air  corridors  are  virtually  unlimited,  but  the  number  of  flights 
that  can  be  assigned  to  routes  is  limited  by  the  number  of  aircraft 
and  crews  available.  For  an  aircraft  constraint,  b  may  represent 
the  total  flying  hours  available,  and  A  the  number  of  flying  hours 
required  for  an  aircraft  to  fly  over  a  c  (i,j).  For  a  crew  constraint, 
b  may  represent  the  number  of  crew  flying  hours  available  and  A„  the 
flying  hours  required  for  aircraft  and  crew  to  fly  over  arc  (i,j). 

The  special  case  of  one  resource  constraint  for  a  single  commodity 
network  flow  has  been  treated  before.  It  is  usually  attacked  by 
assigning  each  arc,  (i,j),  a  cost  of  A^ ,  generating  a  profile  of  flow 
versus  required  cost,  and  choosing  the  point  at  which  cost  is  equal  to  b 
[13-15],  which  would  be  more  efficient  than  the  method  presented  here. 
However,  the  method  of  this  Memorandum  has  the  advantages  that  more 
than  one  resource  constraint  may  be  handled,  multicommodity  networks 
treated,  and  the  resource  constraints  separated  from  the  cost  or 
objective  function. 

Multimode  Networks 

A  special  case  of  resource  constraints  that  deserves  special 
attention  is  the  situation  where  several  modes  of  transport  are  in¬ 
volved  and  one  or  more  vehicle  types  is  in  limited  supply.  No  more 
than  one  vehicle  type  is  allowed  to  carry  units  of  flow  on  any  one 
arc.  This  is  really  no  restriction  because  a  physical  arc  using  more 
than  one  type  of  vehicle  may  be  broken  up  into  several  arcs,  each 
of  which  uses  a  single  type.  The  number  of  vehicle  units  (i.e.,  hours. 
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days,  etc.)  of  a  given  type  consumed  by  an  arc  using  that  vehicle  is 
a  linear  function  of  its  flow.  Thus,  if  a  unit  of  flow  on  arc  (i,j) 

r 

consumed  units  of  vehicle  type  r,  and  the  number  of  vehicle  r 

units  is  b  ,  then  one  would  have  [  A, .  y.  .  s  b  .  where  the  summation 
r’  ij  'ij  r’ 

is  taken  over  all  arcs  using  vehicle  type  r. 

These  types  of  constraints  arise  in  interdiction  problems  be¬ 
cause  the  user  of  a  lines-of-communication  network  may  often  have 
several  modes  of  transport  available,  but  is  often  limited  by  a  shortage 
of  the  different  types  of  vehicles  rather  than  by  actual  network  through¬ 
put  capacity  [11]. 
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The  full  representation  of  program  (2)  would  require  a  complete 
enumeration  of  all  chains  for  the  Q  commodities.  The  number  of  chains, 
even  for  relatively  small  networks,  is  large  enough  to  make  such  enu¬ 
meration  impractical;  also,  an  extension  of  Ford  and  Fulkerson  [1] 
makes  it  unnecessary.  Specifically,  at  least  one  optimal  solu¬ 
tion  to  (2)  is  a  basic  solution,  for  which  the  flow  on  no  more  than  m 
chains  is  nonzero.  Furthermore,  the  simplex  method  of  solution  pro¬ 
ceeds  from  one  basic  solution  to  another,  each  time  dropping  one  column 
and  adding  another  to  the  basis,  until  optimality  is  reached.  Thus, 
at  any  point  of  the  calculation,  it  is  only  necessary  to  maintain  the 
columns  of  the  current  basis,  and  to  be  able  to  generate  the  incoming 
column  for  the  next  basis.  As  will  be  shown  later,  the  generation  of 
the  new  column  can  be  done  by  attaching  appropriate  costs  or  lengths 
to  the  arcs,  and  finding  shortest  chains  for  the  Q  commodities. 

Let  B  be  the  current  basis  and  n  the  corresponding  simplex  multi¬ 
pliers,  Then  multiplying  (2)  by  B  \  one  obtains  for  the  cost  coeffi¬ 
cient  of  x.,  the  quantity 


Ct  + 

-rr  for  t  >  n 
r 

where  rrr  is  the  simplex  multiplier  for  row  r.  From  (3)  and  (4),  one 


m 

/  tt  3  for  t  s  n 
/  y  r  rt 

r-1 


obtains 
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(5.1) 


where 

(5.2) 


Thus,  if  one  were  to  assign  each  arc,  (i,j),  a  length  of  jg ^ ^ ,  cfc  would 
merely  be  the  total  length  of  chain  Ct<  The  current  basis  would  be 
optimal  if  all  tt^  2  0  and  the  length  of  all  chains  were  nonnegative. 
Otherwise,  one  may  perform  a  pivot  operation  on  a  column,  r.  r,  where 
t^  <  0,  or  on  a  column,  j  s  n,  where  the  length  of  is  strictly 
negative,  to  obtain  the  next  basis.  The  search  for  negative  chains  may 
be  accomplished  by  finding  the  shortest  chain  for  each  commodity.  Thus, 
a  generalized  minimum  cost  multicommodity  flow  may  be  found  by  the  fol¬ 
lowing  algorithm, 

1.  Use  as  an  initial  starting  basic  solution,  xn+r  =  b^ 
all  other  xt  =  0, 

2.  Let  tt  be  the  vector  of  the  simplex  multipliers  for  the 
current  basis. 


3.  Let  tt  ■  min  tt  •  If  n  <0,  introduce  x  ,  into  the  basis 

q  r  r  q  n+q 

by  pivoting  and  return  to  2.  Otherwise  continue, 

4.  Assign  each  arc  (i,j)  a  length  equal  to  that  in  expression 

(S.2).  For  each  commodity,  q,  find  the  shortest  chain  C.,  N 


between  its  source  and  sink,  and  denote  its  length  by 
ini.n 

5.  Let  i-  m  i.,  If  C.,  .  i  0,  terminate  as  the  program 

q  q  j(q)  j(q) 

is  optimal.  Otherwise,  introduce  xj^qj  into  the  basis  by  pi¬ 
voting,  and  then  go  to  Step  2, 
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There  are  several  algorithms  available  for  finding  the  shortest 
chain  in  a  network,  and  many  of  these  have  been  examined  and  compared 
by  Dreyfus  [16].  The  first,  and  one  of  the  most  efficient,  is  one  by 
Dykstra  [17].  Hu  [18  ]  has  some  interesting  procedures  for  finding 
shortest  chains  between  all  pairs  of  nodes,  which  may  be  quite  useful 
when  the  number  of  commodities  is  large.  Many  of  these  procedures 
require  that  the  network  contain  no  cycles  of  negative  length.  This 
condition  is  assured  in  Step  4;  because  the  cost  portion  (k  's)  of  all 
cycles  are  nonnegative  by  assumption,  all  a_  ^  0,  and  Step  3  assures 

all  tt  ^  0  before  Step  4  is  entered, 

r 

There  are  several  possible  minor  modifications  to  the  algorithm 
that  may  prove  desirable  and  are  easy  to  make.  One  possibility  is  to 
avoid  finding  shortest  chains  for  all  commodities  by  terminating  the 
search  as  soon  as  a  negative  chain  is  found.  If  this  is  done,  one 
may  wish  to  examine  the  commodities  in  cyclical  order.  Thus,  if 
commodity  q  is  the  last  one  examined  in  one  iteration,  commodity  (q+1) 
mod  Q  is  the  first  one  examined  for  the  next  iteration.  Another  pos¬ 
sibility  is  to  add  the  shortest  chain  columns  for  all  commodities  to 
the  basic  columns,  optimize  the  resulting  program,  and  then  drop  all 
columns  that  do  not  appear  in  the  new  optimal  basis. 
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V,  AN  ALTERNATIVE  METHOD  OF  SOLUTION 

The  algorithm  presented  in  the  last  section  finds  a  generalized 
multicommodity  flow  by  solving  an  arc  subset-chain  linear  program  using 
simplex  multipliers  to  generate  the  new  column  to  enter  the  basis.  The 
column  generated  is  that  with  the  most  negative  cost  coefficient.  Of 
course,  one  could  choose  any  variable  with  a  negative  cost  coefficient 
to  enter  the  basis.  Kuhn  and  Quandt  [19]  have  performed  computational 
experiments  on  several  different  column  generating  schemes  for  general 
linear  programs;  the  most  efficient  ones  proved  to  be  those  which  chose 
the  column  with  the  greatest  negative  cost  coefficient  after  the  columns 
had  been  normalized.  This  normalization  could  be  accomplished  by  div¬ 
iding  the  column  elements  by  their  absolute  sums,  by  the  sum  of  their 
positive  elements,  by  the  sum  of  the  squares  of  the  elements,  or  other 
similar  quantities.  While  many  of  these  schemes  would  require  one  to 
enumerate  all  columns,  modifications  of  them  allow  one  to  avoid  this. 

The  criterion  proposed  here  is  a  modification  of  that  by  which  normal¬ 
ization  is  performed  by  dividing  the  cost  coefficient  by  the  sum  of  the 
absolute  values  of  all  column  elements. 

Specifically,  let  B  be  the  current  basis  and  n  ♦‘he  corresponding 
vector  of  simplex  multipliers,  which  is  also  the  last  row  of  B  Fur¬ 
thermore,  designate  the  element  in  row  r  and  column  t  of  B  *  bv  B  *  . 

J  rt 

Then,  for  the  individual  matrix  elements,  a  ,  one  has 

rt 

a  =  N  ' B" 1  a 
rt  /  rk  kt 

k 


(6.1) 
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i 


(6-2)  E1..  ■£(*.,££) 

r  r  k 

-£[  s  (ZX)  ] 

r  (i,j)eCt  k 

Thus,  i  one  were  to  assign  each  arc,  a  length  of 


the  sum  of  the  elements  of  the  column  representing  chain  j  would 
merely  be  the  sum  of  the  chain's  arc  lengths.  Instead  of  normalizing 
the  cost  coefficient  by  dividing  it  by  the  sum  of  absolute  values  of 
its  elements  (i.e.  terms  on  the  left-hand  side  of  (6.2)),  as  is  often 
done,  normalizing  the  cost  coefficient  in  this  paper  will  be  done  by 
dividing  by  the  sum  of  the  absolute  value  of  the  arc  lengths  in  ex¬ 
pression  (7).  Thus  letting 


(8.1) 


XKj  ZX 


(8.2) 


"t-  E  2u 

(i,j)eCt 


the  next  column  or  variable  to  introduce  into  the  basis,  if  the  current 
basis  is  not  optimal,  is  column  s,  where 


c 

s 


N 


s 


min 


(9) 


"71 
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A  test  to  find  whether  a  particular  chain,  s,  for  which  c  <0  satis- 

8 

fics  (9)  is  suggested  by  the  following  theorem. 


»  I  c 

Theorem  1:  Suppose  c  *  0  and  let  K  ■ _ s 

9  _  I  Ks 

satisfies  (9)  if  and  only  if  e  +  KN  ■  0  ■  min 

s  s  ct 

if  this  condition  is  not  met,  any  chain  with  c^ 
and 


I  .  Then  chain  s 

ct  +  KNt#  Furthermore 

+  KNt  <  0  satisfies  cf 


<r  0 


Proofs  Since  K  ■  Cs/N  ,  c  +  KH  “0. 

/  8  S  S 

all  chains  with  c£  *  0  satisfy  c£  +  KNt  >  0. 

Ihen  «t/»t  *  - Kor  s+ra,t 

t/»i  <  %/».  • The"  <  -  k 


ct  <  0  and  c 


Also,  because  all  f  .  >  0, 

ij 

Suppose  cfc  <  0  and 
;>  0.  Similarly,  suppose 

or  c  +  KUt  <  0.  QED 


Noting  that  if  one  assigns  arc  ( i ,  j)  a  length  of  j l  +  K?.^, 

the  network  contains  no  negative  cycles,  and  the  length  of  chain  Ct 
is  ct  +  KNt,  it  follows  that  Steps  4  and  5  of  the  algorithms,  which 
find  the  column  to  enter  the  basis  and  determine  when  optimality  is 
reached,  may  be  modified  as  follows. 

4.  For  each  commodity  q,  do  the  following: 

(a)  Assign  each  arc  (i,j)  a  length  as  defined  in  ex¬ 
pression  (5.2),  and  find  the  shortest  chain.  If  its  length  is  non- 
negative,  designate  it  as  and  go  to  5.  Otherwise  designate 

it  as  C  . 
m 


t 


I 


u 
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(b)  Let  K  ■  c  /N  where  I  c  and  N  are  as  defined  in 
m|/  m  |  m  m 

(5.1)  and  (6.2).  Then  assign  each  arc  (i,j)  a  length  of 
and  find  the  shortest  chain. 


(c)  If  the  length  of  the  shortest  chain  found  is  nonnegative, 

let  t (q)  =  m,  and  go  to  5.  Otherwise,  let  the  shortest  chain  found 

replace  C  , and  go  back  to  4(b). 
m 

5.  If  min  c  ,  .  *  0,  terminate,  as  the  current  basic  solution 
q  t(q) 

is  optimal.  Otherwise,  let  q  be  such  that  cfc  (q)/Nt(q)  “  T  Ct(q)/Nt(q)» 


introduce  into  basis,  and  go  back  to  Step  2. 

As  was  the  case  with  the  other  criteria,  the  procedure  may  be 


modified  to  select  a  chain  from  the  first  commodity  found  to  have  a 


negative  length  chain  and  to  consider  commodities  in  cyclical  order. 
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VI.  NODE-ARC  FORMULATION  AND  DECOMPOSITION 

The  generali. jd  multicommodity  flow  problem  may  be  formulated 
in  terms  of  its  node-arc  incidence  matrix,  as  well  as  its  arc  subset- 
chain  incidence  matrix.  Specifically,  let 

"'ij  =  flow  of  commodity  q  on  arc  >i,j) 
r^  *  the  total  flow  on  commodity  q  chains 
Then  the  multicommodity  flow  problem  is  defined  by  the  following  con¬ 
straints. 


Find  min  z,  y^  2.  0, 


r^  ^  0  such  that 


(10.1) 

2-  EEkuyu 

I  <i,j) 

(10.2) 

r'1* 

q  uTj) 

(10.3) 

EVE’!*  • 

j  j 

r  i  =  s 

q  q 

0  i  J  s  ,  t 

f  q  q 

■  -r  i  =  t 

q  q 

for  q  =  1 ,  ....  Q 

The  transportation  structure  of  (10.3)  allows  one  to  attack  this 
program  efficiently  by  the  Dantzig-Wolfe  decomposition  algorithm  [20], 
with  (10,2)  serving  as  the  master  program  and  the  Q  sets  of  equations 
of  (10.3)  forming  Q  subprograms. 

Note  that  all  constant  terms  in  the  Q  subprograms  are  zero,  and 
consequently  none  of  these  programs  have  extreme  point  solutions  other 
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than  aLL  y*? .  =  0.  Thus,  if  one  lets  (y*?!),  .... 

lj  ij  ij 

vectors  spanning  the  solution  space  for  subprogram  q,  (10) 
formulated  as  follows  when  decomposition  is  applied. 


be  a  set  of 
is  re  - 


(11.) 


Find  min  z,  x  >  0,  s  ^  0,  v  ^  0  such  that 
qt  r  r 


z  =  >  c  x 

qt  qt 


q.t 


(11.2) 


q,t 


P  x  +v  =b  r  *  1 ,  . . .  ,  m 

qt  qt  r  r 


where 


(12.1) 

Pr  =  VAr. 

qt 
y?. 
1 J 

(i.  j) 

(12.2) 

c  =  y  k. . 
qt  Z— *  tJ 

-qt 

y  •  • 
ij 

(i>  j) 

Note,  however,  that  any  feasible  solution  to  subprogram  q  (10.3) 

may  be  expressed  as  a  nornegative  linear  combination  of  commodity  q 

chain  flows  and  that  solutions  corresponding  to  distinct  chain  flows 

are  linearly  independent.  Thus,  the  (y^)  may  be  defined  more  speci- 

•O 

fically  as  the  N(q)  solutions,  corresponding  to  one  unit  of  flow  on  a 


commodity  q  chain  (i.e.  r^  =  1).  When  this  is  done,  program  (11)  is 

identical  to  program  (2).  Specifically,  the  c  and  P  of  (11)  are 

equal  to  the  c.  and  a,,  of  (2)  respectively,  and  the  x  .  and  s  of  (11) 
j  ij  r  qi  r  v  ' 


* 

Strictly  speaking,  such  a  solution  may  also  contain  cycles. 
However,  as  cycles  are  nonnegative  by  assumption,  they  cannot  appear 
in  an  optimal  solution,  and  consequently  need  not  be  considered. 


* 
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are  identical  to  the  x^  of  (2).  Furthermore,  when  the  Dantzig-Wolfe 
decomposition  is  applied,  not  all  columns  of  (11)  are  enumerated 
during  the  course  of  the  calculation.  Instead,  only  those  correspon¬ 
ding  to  the  current  basis  are  maintained.  Letting  nr  be  the  simplex 
multiplier  corresponding  to  row  r  for  the  current  basis,  the  next 
column  to  enter  the  basis  (i.e.  that  with  the  lowest  cost  coefficient 

cqt  'XX  Pqt  ■  XXlj  +  XXr  *U  1  >  if  a11  "r  *  0  ls  £ou‘'d 

r  (i,  j)  r 

by  solving  the  following  Q  programs.  If  z-  ■  min  (min  z  )  <  0,  then 

q  q  q 

the  new  column  is  generated  by  applying  the  transformation  (12)  to 
the  solution  for  program  q. 

Find  min  z  .  y? .  ^0  such  that 

q  ij 


(13.1) 

(13.2) 


However,  if  arcs  are  assigned  lengths  equal  to  the  I  of  (5,2),  an 

optimal  solution  to  (13)  may  be  obtained  by  finding  a  shortest  path 

from  8^  to  tq  and  setting  y^  -  1  if  (i,j)  is  on  this  path  and  y^  -  0 

otherwise.  Of  course,  if  some  tt  <0,  then  v  is  the  variable  to  enter 

r  r 


the  basis.  The  procedure  terminates  when  all  tt  i  0  and  all  min  z  a  0. 

r  q 


The  above  procedure,  however,  is  identical  to  that  given  in 


Sec.  III.  Thus  applying  the  Dantzig-Wolfe  decomposition  algorithm 
to  the  node-arc  formulation  is  identical  to  using  the  arc  subset-chain 
procedure  given  earlier. 
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VII.  FLOW  MAXIMIZATION  AND  FEASIBILITY 

Two  multicommodity  flow  problems  that  often  arise  in  practice  are 
that  of  maximizing  a  linear  combination  of  the  commodity  flows,  and 
that  of  finding  a  feasible  routing  which  will  meet  required  flows. 

This  section  shows  how  these  problems  may  be  formulated  in  terms  of 
finding  generalized  minimum  cost  flows. 

Maximizing  a  Linear  Combination  of  the  Commodity  Flows 

To  maximize  a  linear  combination  of  commodity  flows,  say  a^P^  + 
^2^2  +  •••  +  3qPq  where  P^  is  the  total  flow  of  commodity  i,  one  need 
only  attach  an  artificial  node,  s^  for  each  commodity,  and  an  artifi¬ 
cial  arc  directed  from  s^  to  node  i's  true  source,  and  assign  it  a 
cost  of  -a^,.  It  can  easily  be  seen  that  if  all  other  arcs  have  zero 
costs  and  essentially  infinite  upper  flow  bounds  on  each,  the  cost 
minimization  problem  is  the  same  as  the  linear  combined  flow  maximiza¬ 
tion  problem.  If  other  arc  costs  are  not  zero,  then  care  must  be  taken 
to  insure  that  the  negative  of  any  artificial  arc  cost  is  greater  than 
the  length  of  any  chain.  This  may  be  accomplished  by  multiplying  all 
the  a^'s  by  a  sufficiently  large  positive  constant.  It  should  be  noted 
that  none  of  the  artificial  arcs  belong  to  cycles,  and  therefore  their 
addition  will  not  cause  the  restriction  of  no  negative  cycles  to  be 
violated.  It  should  also  be  noted  that  since  these  arcs  really  have 
no  capacity  constraints,  their  addition  adds  no  arc  subset  constraints, 
and  hence  no  constraints  to  the  linear  program.  Essentially  their  only 
effect  is  to  terminate  the  program  in  Step  5,  with  each  l . ,  .  ^  -a  in- 

J  (q)  *  °' 


stead  of  t 
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Meeting  Required  Flows 

The  problem  of  meeting  required  commodity  flows,  P  , 

may  be  accomplished  by  adding  the  same  artificial  arcs  and  nodes  as 
when  maximizing  a  linear  combination  of  commodity  flows.  However, 
the  artificial  arc  connected  to  must  be  assigned  an  upper  flow 
bound  or  capacity  of  P^  and  a  cost  sufficiently  negative  to  insure 
that  all  chains  for  commodity  i  have  negative  length.  For  this  prob¬ 
lem,  each  of  the  artificial  arcs  does  add  a  linear  constraint  to  the 
program. 
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(i.j) 


k.  . 
ij 
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q 

t 

q 
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b 
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LIST  OF  SYMBOLS 

the  arc  directed  from  node  i  to  node  i. 
erst  of  a  unit  of  flow  on  arc  (i,j). 
source  node  for  commodity  q. 
sink  node  for  commodity  q. 
the  number  of  commodities. 

upper  bound  or  capacity  for  linear  flow  constraint  r. 

contribution  of  a  unit  of  flow  on  arc  (i,j)  to  flow 
constraint  r. 

total  flow  on  arc  (i,j). 

flow  of  commodity  q  on  arc  (i,j). 

r 

subset  of  arcs  for  which  A. .  >0. 

ij 

network  commodity  chain  t. 
sum  of  arc  costs  on  chain  t, 

contribution  of  a  unit  of  flow  on  chain  t  to  flow  constraint  r. 
current  basis  for  the  linear  program, 
simplex  multiplier  for  row  r. 

length  of  arc  (i,j)  as  determined  by  arc  costs,  constraint 
coefficients,  and  simplex  multipliers. 

total  chain  t  cost  relative  to  the  basis  B. 

the  coefficient  matrix  element  in  row  r  column  (chain)  t 
with  respect  to  the  basis  B 

the  row  r,  column  t  element  of  B 

length  of  arc  (i,j)  as  determined  by  the  coefficient  matrix 
(not  cost  row)  with  respect  to  basis  B. 


total  flow  of  commodity  q 
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product  of  joint  constraint  row  r  and  a  feasible  solution 
to  the  node-arc  incidence  program  for  commodity  q. 

product  of  the  cost  row  of  the  multicommodity  program  and 
a  feasible  solution  to  the  node-arc  incidence  program 
for  commodity  q. 
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